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ABSTRACT 


A  comprehensive  field  rip-current  experiment  (RCEX)  was  conducted  from  14 
April  to  18  May  2007  in  Sand  City,  California,  on  Monterey  Bay.  Lagrangian 
observations  were  made  with  inexpensive  ($150),  handheld,  Differential  Global 
Positioning  Systems  (DGPS)  mounted  on  surf  zone  drifters.  The  inexpensive  DGPS 
requires  post-processing  to  achieve  accuracy  0(0. 4m)  and  0(0.0 1mA).  Thirty  drifters 
were  constructed  and  deployed  in  well-developed,  rip-currents  to  map  the  circulation 
patterns  for  the  first  time  in  the  field  at  a  high  spatial  resolution.  Drifter  observations 
obtained  during  three-hour  periods  on  seven  different  days  under  varying  wave  and  tidal 
conditions  describe  eddies  with  a  rotational  period  of  4.7mm,  confined  to  the  surf-zone 
and  coupled  to  the  rip  morphology.  On  average,  three  drifters  per  hour  exited  the  surf- 
zone.  Dependent  upon  wave  conditions,  one  or  two  eddies  existed  between  90m-spaced 
rip-channels,  creating  a  seaward  flow  in  the  channels  and  shoreward  flow  over  the  shoals. 
Cross-shore  volumetric  flow  rates  for  an  alongshore  transect  through  the  eddy  centers 
balance  to  a  difference  of  less  than  10%  of  the  gross  flow  discharge.  Velocity 
measurements  obtained  from  drifter  data  are  evaluated  with  velocities  obtained  from 
stationary,  in-situ  instruments. 
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I.  INTRODUCTION 


A.  MOTIVATION 

The  ocean  has  potential  to  both  give  life  and  take  it  away,  especially  near  shore. 
According  to  the  United  States  Lifesaving  Association  (USLA),  100  deaths  occur 
annually  due  to  rip  currents,  and  80%  of  all  beach  lifeguard  rescues  are  rip  current  related 
(USLA,  2007).  There  is  a  great  amount  of  interest  in  surf  zone  circulation  with  a  wide 
variety  of  applications,  including  sediment  and  pollution  transport,  shoreline  evolution, 
and  beach  and  swimmer  safety.  Furthermore,  research  focused  on  rip  currents  is 
important  in  predictive  surf  zone  modeling  that  ultimately  aids  military  planners  in  mine¬ 
hunting  and  SEAL  delivery. 

Rip  currents  are  typically  “shore-normal,  narrow,  seaward-flowing  currents”  that 
extend  from  the  shoreline  to  the  edge  of  the  surf  zone  with  relatively  high  speeds 
(MacMahan  et  al.,  2005).  Since  the  1930s,  field  experiments  have  collectively  improved 
our  understanding  of  rip  current  dynamics  (Shepard  et  al.,  1941;  Shepard  and  Inman, 
1950;  Sonu,  1972;  Brander,  1999;  Brander  and  Short,  2000;  Schmidt  et  al.,  2005; 
MacMahan  et  al.,  2005).  Owing  to  the  cost,  difficulty,  and  danger  of  installing  multiple 
instruments  throughout  a  rip  current  system  many  hypotheses  remain  unanswered. 
Lagrangian  observations  have,  for  the  most  part,  qualitatively  complemented  Eulerian 
measurements,  and  have  the  potential  to  resolve  these  rip  current  hypotheses.  Until  now, 
the  cost  of  accurate  drifter  measurements  has  limited  the  numbers  deployed,  and 
therefore  our  ability  to  quantitatively  understand  rip  current  systems.  The  goal  of  this 
thesis  is  to  transform  previously  qualitative  Lagrangian  descriptions  of  rip  currents  into  a 
more  quantitative  description,  using  a  fleet  of  inexpensive,  accurate,  self-logging  drifters 
recently  deployed  in  the  field. 

B.  RIP  CURRENT  DEFINITION 

A  rip  current  can  essentially  be  thought  of  as  a  localized  river  in  the  ocean, 
carrying  water  from  the  shoreline  out  to  sea.  Differences  in  bathymetry  lead  to  variations 
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in  wave  breaking  along  shore  that  cause  a  higher  mean  water  level  over  shoals  than  over 
rip  channels.  The  resulting  circulation  pattern  (Figure  1)  is  offshore-directed  flow  in  the 
rip  channel  and  onshore  flow  over  the  shoals  to  either  side. 


Figure  1.  Conceptual  representation  of  a  rip  current  circulation  system.  Shaded 
regions  represent  shoals  to  either  side  of  a  rip  channel.  Wavy,  dotted  line 
represents  line  of  wave  breaking.  Onshore  flow  occurs  over  the  shoals.  Onshore 
flow  is  funneled  through  feeder  channels  into  the  rip  channel  and  seaward  as  part 
of  the  rip  current.  [After:  MacMahan  et  ah,  2005]. 

A  conceptual  drawing  illustrates  the  primary  aspects  of  a  rip  current  cell  (Figure 
1).  The  dotted  line  represents  wave  breaking  -  0(1 00m)  from  the  shoreline  in  the  recent 
rip  current  experiment  (RCEX),  but  this  is  dependent  on  wave  and  tide  conditions.  When 
most  beach-goers  describe  rip  currents,  they  think  of  the  rip  neck,  the  narrowest  section 
of  the  rip  characterized  by  strong  seaward  flow.  The  rip  head  is  the  seaward  tip  of  the 
rip,  where  rip  currents  diffuse  and  lower  velocities  are  present.  The  flow  in  the  rip  head 
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is  less  constrained  by  bathymetry  of  the  rip  channel  and  spreads  laterally.  Returning 
onshore  flow  is  funneled  into  the  rip  through  alongshore-directed  feeder  currents,  which 
occur  near  the  shoreline. 

In  1941,  Shepard  et  ah,  showed  the  intensity  of  rip  currents  increases  with 
increasing  wave  height  and  observed  a  similar  trend  in  their  seaward  extent.  From  this 
they  concluded  that  rip  currents  result  from  wave  translational  motion  piling  up  along  the 
shoreline,  which  must  return  offshore  to  satisfy  continuity.  Some  of  the  returning  water 
is  concentrated  in  zones,  which  develop  in  persistent  outward  currents,  called  rip 
currents.  This  hypothesis  and  several  others  have  been  proposed  to  explain  the  formation 
of  rip  currents. 

Other  potential  mechanisms  for  rip  current  generation  have  been  separated  into 
two  main  categories:  wave  interaction  and  structural  (boundary)  interaction  (Dalrymple, 
1978).  The  wave  interaction  theories  suggest  that  rip  currents  are  caused  and 
strengthened  by  variations  in  mean  water  level  along  the  shore  leading  to  hydrostatic 
differences.  Such  theories  include:  interaction  between  incident  edge  waves  (traveling 
along  the  shore)  (Bowen,  1969;  Bowen  and  Inman,  1969);  intersecting  wave  trains 
(Dalrymple,  1975);  or  wave-current  interaction  (LeBlond  and  Tang,  1974;  Dalrymple  and 
Lozano,  1978).  Boundary  interaction  theories  suggest  that  rip  currents  are  the  result  of 
circulation  patterns  created  by  wave  field  interaction  with  topographical  variations,  such 
as  bathymetry  (Bowen,  1969;  Noda,  1974;  Dalrymple  et  ah,  1976)  or  artificial  structures 
Liu  and  Mei,  1976;  Mei  and  Angelides,  1977)  -  e.g.  jetties.  Problems  with  these  theories 
were  that  the  fonner  under-predicted  surf  zone  energy  dissipation,  while  the  latter  failed 
to  accurately  predict  the  wave  field  over  longshore  sand  bars  (Dalrymple,  1978). 

C.  HISTORICAL  CONTEXT 

Surf  zone  circulation  and  dynamics  have  been  studied  for  several  decades,  and 
interest  in  rip  currents  can  be  traced  back  to  the  mid- 193 Os  when  a  scientific  debate  about 
the  existence  of  an  “undertow”  led  Shepard  (1936)  to  coin  the  tenn  “rip  current”  as  a 
descriptive  name  for  the  phenomenon  that  seemed  to  vex  lifeguards.  Subsequent 


3 


research  has  shown  that  undertow  does  exist  (Hansen  and  Svendsen,  1984),  but  occurs  on 
two-dimensional  beaches  rather  than  three-dimensional  beaches  where  rip  currents  are 
the  norm. 

Early  experiments  were  primarily  qualitative,  visually  tracking  dye,  wave 
breaking  patterns,  kelp,  and  other  floating  objects  (Shepard,  1936;  Shepard  et  ah,  1941; 
Shepard  and  Inman,  1950).  Rip  currents  in  these  studies  were  characterized  by  velocities 
from  9cm/s-22cm/s  (Shepard  et  ah,  1941);  widths  from  15-150/u  and  offshore  extents 
from  100-800/n  (Shepard,  1936;  Shepard  et  ah,  1941);  and  spacing  O(180m).  These  early 
studies  generalized  nearshore  circulation  as  net  shoreward  flow  in  wide  lanes  between 
seaward  flowing  rip  currents  (Shepard  and  Inman,  1950)  and  permitted  qualitative 
representations  of  nearshore  flow  (Figures  2  and  3). 


Figure  2.  Sketch  of  a  rip  current  system  based  on  qualitative  estimates  of  flow  from 
early  field  experiments.  [From:  Shepard  et  ah,  1941]. 
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Figure  3.  Detailed  qualitative  circulation  patterns  in  the  nearshore  as  detennined 

from  field  experiments  in  the  1940s.  [From:  Shepard  and  Inman,  1950]. 

Later  experiments  compared  visually  derived  measurements  with  measurements 
from  static  instruments  in  the  rip  channels  and  on  the  shoals  (Sonu,  1972;  Brander  and 
Short,  2000).  These  experiments  provided  more  quantitative  pictures  of  surf  zone 
circulation  (Figures  4  and  5).  Circulation  was  characterized  by  mean  rip  current 
velocities  0(60-70cmA)  (Sonu,  1972;  Brander  and  Short  2000),  with  a  maximum  as  high 
as  2 m/s  (Sonu,  1972),  and  onshore  flow  between  rip  currents  from  20-30 cm/s  (Sonu, 
1972). 
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Figure  4.  Flow  field  derived  from  measurements  in  the  early  1970s.  The  seaward 

flowing  rip  current  can  be  seen  where  the  streamlines  are  closest  together  (center), 
flanked  on  either  side  by  shoreward  flow  with  more  widely  spaced  streamlines. 
[From:  Sonu,  1972]. 
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Figure  5.  Lagrangian  measurements  of  a  rip  channel,  using  floats  tracked  by  two 
shore-based  theodolites.  Lack  of  circulation  likely  due  to  rip  channel  angled 
obliquely  from  shore.  [From:  Brander  and  Short,  2000], 
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More  recently,  Global  Positioning  System  (GPS)-tracked  drifters  have  further 
improved  our  understanding  of  the  surf  zone  (Schmidt  et  al.,  2005),  to  show  rip  current 
intensification  as  tide  level  decreases.  Unfortunately,  owing  to  the  high  cost,  0($6K),  of 
these  systems,  only  10  drifters  were  used,  and  drifters  were  only  released  in  pairs, 
providing  a  limited  physical  picture  of  the  flow  structure  (Figure  6). 
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Figure  6.  Quantitative  measurements  of  rip  channels,  obtained  using  10  GPS- 
tracked  surf  zone  drifters  deployed  in  pairs  during  five-  to  six-hour  periods. 

[From:  Schmidt  et  ah,  2005]. 

Recent  laboratory  experiments  using  Acoustic  Doppler  Velocimeters  (ADVs)  to 
measure  current  patterns  on  a  barred  beach  with  rip  channels  suggest  that  bars  and  rip 
channels  dominate  nearshore  dynamics  (Haller  et  ah,  2002).  These  results  provide  good 
physical  detail  and  show  distinct  eddies  with  offshore  flow  over  rip  channels  and  onshore 
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flow  over  adjacent  shoals  (Figure  7).  However,  the  laboratory  conditions  used  in  this 
experiment  do  not  accurately  reflect  field  conditions,  since  the  waves  were 
monochromatic  (Haller  et  ah,  2002).  Laboratory  experiments  often  use  a  smaller  scale 
than  field  conditions,  and  are  thus  limited  in  their  ability  to  reproduce  realistic  conditions. 
Additionally,  laboratory  results  are  often  difficult  to  reproduce  under  realistic  field 
conditions,  owing  to  the  cost  and  difficulty  of  placing  such  a  large  array  in  the  field. 


Figure  7.  Mean  flow-field  from  a  laboratory  tank  using  monochromatic  waves  on  a 
barred  beach,  obtained  by  moving  two  instruments  from  point-to-point. 

Circulatory  patterns  are  observed  with  onshore  flow  over  the  bars  and  offshore 
flow  through  the  channel  between  adjacent  bars.  [From:  Haller  et  ah,  2002], 

While  the  numerical  results  focusing  on  rip  currents  have  certainly  increased  our 
understanding  of  surf  zone  circulation,  the  results  are  largely  evaluated  at  a  limited 
number  of  fixed  instrument  arrays  that  have  not  captured  the  full  development  of  rip 
current  circulation  cells.  Until  recently,  the  scarcity  of  Lagrangian  data  has  meant 
Eulerian  measurements  have  been  the  primary  data  available  for  predictive  surf  zone 
modeling  (Johnson  and  Pattiaratchi,  2004). 

Furthermore,  owing  to  the  high  cost  of  GPS-tracked  surf  zone  drifters,  most 
Lagrangian  observations  have  taken  place  in  the  laboratory  under  ideal  bathymetric 
conditions  (Kennedy  and  Thomas,  2004),  and  those  that  have  taken  place  in  the  field 
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have  been  limited  in  size  (Schmidt  et  al.,  2005).  Recently,  an  inexpensive  ( — $  1 50)  hand¬ 
held  GPS  device  was  tested  for  use  in  oceanographic  drifters  (MacMahan  and  Brown, 
2007).  The  positive  results  of  these  tests  permitted  construction  of  a  larger  fleet  of 
drifters  (30  drifters)  for  the  2007  rip  current  experiment  (RCEX)  at  a  significantly  lower 
cost,  without  sacrificing  precision. 
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II.  EXPERIMENTAL  DESCRIPTION 


A.  SUMMARY 

A  field  rip-current  experiment  (RCEX)  was  conducted  from  14  April  to  18  May 
2007  to  examine  rip-current  circulation  patterns  in  greater  detail.  Eulerian  measurements 
were  obtained  using  a  stationary  array  of  1 1  instruments.  Lagrangian  measurements  were 
obtained  with  30  inexpensive  drifters. 

B.  FIELD  SITE 

1.  Site  Description 


Figure  8.  (a)  Plan  view  of  field  site  used  for  the  Rip  Current  Experiment  (RCEX). 

(b)  Close-up  of  field  site  with  sign  warning  bathers  about  submerged  equipment. 

The  field  site  was  approximately  a  300m  alongshore  by  250m  cross-shore  region 
located  in  Sand  City,  California,  at  the  south  end  of  the  Monterey  Bay  (Figure  8). 
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Bathymetric  profile  along  Y=  115m  for  5/04/2007 


Figure  9.  Cross-shore  bathymetric  profile  over  the  shoal  at  the  site.  Zero  datum  in 
elevation  corresponds  to  mean  sea  level  for  4  May  2007.  Cross-shore  coordinates 
show  meters  from  experimental  origin  rather  than  meters  from  shoreline. 


This  is  the  same  site  used  in  the  2001  RIPEX  experiment  and  is  characterized  by 
shore-normal  waves,  coarse  sand,  and  persistent,  energetic  rip  currents  (MacMahan,  et 
ah,  2005).  The  beach  foreshore  is  steep  (1:9),  with  a  low-tide  terrace  (1:100),  and  an 
offshore  slope  of  1:20  (Figure  9). 

For  consistency,  the  local  coordinate  system  is  also  based  on  RIPEX  coordinates 
(MacMahan,  et  al,  2005).  The  mean  shoreline  (zero  vertical-axis  datum  level)  is  about 
50 m  from  the  origin,  and  parallel  with  the  alongshore  component  of  the  array.  Positive  is 
toward  the  North  and  offshore. 
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2. 


Static  Array 
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Figure  10.  Static  array  and  site  bathymetry.  X-axis  is  alongshore  direction  and 
positive  to  the  left.  Y-axis  is  cross-shore  direction  and  positive  up.  Red  dots 
represent  in  situ  instruments  and  are  plotted  in  the  local  coordinate  system  against 
bathymetry  typical  of  the  field  site.  Mean  shoreline  is  the  border  between  blue 
contours  (submerged  bathymetry)  and  yellow  contours  (exposed  beach). 

The  stationary  array  comprised  of  an  intersecting  alongshore-cross-shore  array 
(Figure  10).  The  alongshore  array  consisted  of  seven  digital  electromagnetic  current 
meters  co-located  with  pressure  sensors  (PUVs)  mounted  on  10-foot  poles  that  were  air- 
jetted  approximately  seven  feet  into  the  seabed.  The  cross-shore  component  of  the  array 
consisted  primarily  of  Acoustic  Doppler  Current  Profilers  (ADCPs)  mounted  on  10-foot 
poles  jetted  into  the  bed  approximately  five  feet,  deployed  down  the  center  of  the  main 
rip  channel.  The  instruments  were  inspected  every  couple  of  days,  cleared  of  sand  and 
kelp,  and  raised  or  lowered  as  necessary,  in  order  to  ensure  continuity  and  consistency  of 
measurements.  These  observations  were  part  of  RCEX,  but  are  not  the  focus  of  this 
thesis. 
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3.  Bathymetric  Surveys 

Three  different  types  of  sub-arid  and  sub-aqueous  surveys  were  perfonned  as  part 
of  the  experiment  to  capture  morphologic  changes.  Large  spatial  extent  surveys  of  the 
beach  face  were  perfonned  using  an  All-Terrain  Vehicle  (ATV)  during  low  tide  on  15 
March  and  9  May  2007.  A  wave-rider  (personal  watercraft,  PWC)  was  used  for 
bathymetric  surveys  of  the  site  (from  inside  the  surf  zone  to  200/h  from  the  shoreline) 
during  high  tide  on  20  March,  21  April,  and  1,  11,  and  18  May  2007.  The  presence  of 
waves  and  bubbles  inside  the  surf  zone  made  it  difficult  to  obtain  measurements  by  PWC 
alone.  Walking  surveys  of  the  beach  face  and  shoals  at  the  field  site  were  also  conducted 
during  low  tide  on  20  March,  19,  22,  and  23  April,  and  2,  8,  11,  and  18  May  2007. 
Additionally,  the  instrument  locations  were  surveyed  on  22  April,  and  1 1  and  1 8  May 
2007.  A  total  of  four  complete  surveys,  and  three  instrument  location  surveys  were 
performed.  The  combined  ATV,  PWC,  and  walking  surveys  overlapped  enough  to 
provide  a  good  picture  of  beach  and  surf  zone  morphology  in  the  area  of  interest. 


Figure  11.  Static  array  plotted  against  bathymetric  surveys  from:  (a)  22  Apr,  (b)  1 
May,  (c)  1 1  May,  and  (d)  18  May  2007.  Spikes  in  bathymetry  (11  and  18  May) 
are  caused  by  noise.  Note  consistent  rip  channel  locations  between  surveys. 
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To  simplify  calculations,  an  assumption  is  made  that  bathymetry  does  not  change 
significantly  between  bathymetric  surveys.  There  is  very  little  change  in  the  bathymetry 
over  the  course  of  several  surveys  (Figure  15),  demonstrating  the  validity  of  this 
assumption.  Consequently,  depth  is  approximated  for  any  point  at  the  field  site  given  the 
most  recent  bathymetric  survey  by  accounting  for  mean  sea  level  measured  during  each 
deployment.  This  process  is  later  described  in  detail. 

C.  DRIFTERS 

Surf  zone  drifters  have  been  used  in  research  for  decades  -  even  volleyballs  and 
kelp  used  in  early  experiments  are  technically  drifters  (Shepard,  et  ah,  1941;  Shepard  and 
Inman,  1950).  Advances  in  technology  pennitted  installation  of  small  GPS  receivers  on 
drifters  that  have  improved  measurement  accuracy  (Johnson,  et  ah,  2003;  Schmidt,  et  ah, 
2003;  Schmidt,  et  ah,  2005).  The  high  cost  of  commercially  available  drifters  using  this 
technology  has  made  it  difficult  to  obtain  more  than  10  for  any  experiment  and 
unpleasant  to  risk  loss  of  any,  resulting  in  limited  detail  (Schmidt,  et  ah,  2005).  Less 
expensive  drifters  O($500)  have  been  built  at  the  expense  of  either  reduced  accuracy  or  a 
time-consuming  process  to  match  the  accuracy  of  those  available  commercially  (Johnson, 
et  ah,  2003). 

The  focus  of  this  research  effort  was  an  inexpensive  drifter  designed  to  use 
handheld,  off-the-shelf  Differential  GPS  (DGPS)  technology.  MacMahan  and  Brown 
(2007)  modified  a  commercially  available  DGPS  logger  O($150)  to  obtain  system 
accuracy  O(0.4m)  in  position  and  0(0.01  m/s)  for  velocity  after  post-processing.  The  low 
cost  (less  than  $250  each  without  labor)  of  the  drifters  built  with  this  technology  allowed 
construction  of  a  fleet  of  30,  for  about  the  same  cost  as  a  single  commercially  available 
drifter.  The  large  number  of  drifters  permitted  data  collection  with  greater  spatial  and 
temporal  resolution  than  achieved  by  previous  studies. 
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Figure  12.  Detailed  diagram  of  drifter  construction. 

Figure  12  is  a  graphical  representation  of  a  drifter.  The  drifter  body  designed  for 
this  experiment  is  a  55 cm  length  of  PVC  pipe  with  an  outer  diameter  of  12cm,  and  is 
surmounted  by  a  75 cm  antenna  pole  of  PVC  pipe  with  an  outer  diameter  of  2.6cm  -  for 
an  overall  height  of  1.3/m.  A  35cm-dia.  disk  of  PVC  is  attached  to  the  bottom  of  the  body 
to  dampen  vertical  motion,  and  27c/M-dia.  octagonal,  stainless  steel  plate  is  attached  to  the 
bottom  of  that  to  add  structural  integrity  and  improve  the  drifter’s  righting  moment. 
Three  fins  (46cm  tall  by  7.5cm  wide)  are  equally  spaced  around  the  body  to  improve  it’s 
flow  response.  The  body  is  filled  with  enough  buckshot  to  give  the  drifter  its  desired 
draft  and  improve  its  righting  moment,  and  expanding  foam  is  added  inside  to  prevent  the 
buckshot  ballast  from  shifting.  A  water-proof  case  is  attached  to  the  bottom  of  the 
antenna  pole  to  house  the  DGPS  logger,  and  an  external  antenna  is  mounted  on  top  of  the 
antenna  pole.  The  total  price  tag  of  $250  per  drifter  (without  labor)  includes:  the  DGPS 
logger,  external  antenna,  waterproof  box,  PVC  to  build  the  drifter,  enough  buckshot  for 
ballast,  and  expanding  foam  to  stabilize  the  buckshot. 
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III.  DATA 


A.  DATA  COLLECTION 

1.  Method 

To  collect  Lagrangian  data,  drifters  were  deployed  on  seven  days:  27  April,  and  4, 
5,  7,  10,  15  and  19  May  2007.  Deployment  periods  lasted  from  approximately  0830  to 
1200  on  each  day.  On  a  given  deployment  day,  five  or  six  people  would  initially  start 
loggers  for  10  to  12  drifters.  Once  the  lights  of  the  loggers  were  flashing,  to  indicate  they 
were  collecting  data,  the  drifters  were  taken  to  the  water,  and  deployed  as  close  to  each 
other  as  possible  without  their  bodies  interfering  with  each  other.  As  soon  as  the  first 
group  of  drifters  were  released  in  the  surf  zone,  another  group  of  10  to  12  would  similarly 
be  started  and  deployed.  Drifters  were  removed  from  the  water  when  they  started  to  drag 
along  shoals  or  when  they  left  the  boundaries  of  the  observational  domain.  Drifters  that 
left  the  surf  zone  offshore  were  picked  up  by  boat  and  turned  off  for  the  day.  Throughout 
the  deployment  period  drifters  were  deployed  in  similar-sized  groups  composed  of  new 
drifters  and/or  those  that  washed  up  on  shore  or  crossed  the  longshore  boundaries.  At 
one  point,  as  many  as  25  drifters  were  simultaneously  circulating  in  the  surf  zone.  At  the 
end  of  the  day,  logged  data  were  uploaded  to  computer  for  post-processing,  and  data  were 
cleared  from  the  loggers  for  the  next  deployment. 

Data  from  the  stationary  array  were  used  for  comparison  to  evaluate  flow. 

2.  Sources  of  Error 

There  were  several  sources  of  error  introduced  to  the  data  during  collection. 
Generally  the  easiest  to  identify  was  human  influence  on  the  drifters.  Since  the  logging 
devices  were  turned  on  before  the  drifters  were  taken  to  the  water,  there  were  a  few 
minutes  of  positional  data  at  the  beginning  of  each  set  that  had  nothing  to  do  with  the 
flow  field.  The  same  is  true  whenever  drifters  were  removed  from  the  water. 
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Additionally,  drifters  would  occasionally  become  trapped  on  shoals  or  shoreline  close 
enough  to  a  rip  current  that  it  was  given  a  push  to  set  it  back  in  the  current.  Times  of 
occurrence  were  recorded  by  a  note -taker. 

Another  source  of  error  was  drifter  entanglement  in  kelp.  Floating  kelp  bodies 
attached  to  a  drifter  would  unfavorably  bias  drifter  speeds  and  had  to  be  removed  as  soon 
as  possible.  Since  bodies  of  floating  kelp  were  fairly  large,  it  was  generally  easy  to 
identify  when  they  entangled  a  drifter.  When  drifter  entanglement  was  noticed,  a 
swimmer  would  remove  the  kelp  from  drifter  as  soon  as  practicable,  and  the 
entanglement  and  disentanglement  times  were  recorded  by  a  note-taker. 

Perhaps  the  most  difficult  sources  of  error  to  identify  were  breaking  waves. 
While  the  drifters  are  designed  to  prevent  excessive  vertical  oscillations  due  to  waves 
(fins  and  a  vertical  damper),  breaking  waves  affected  the  drifters  by  occasionally 
submerging  their  antennas.  Drifters  were  designed  to  pass  under  oncoming  waves  (to 
prevent  surfing)  and  their  signal  was  periodically  interrupted  or  degraded  when  this 
occurred.  Also,  during  the  first  deployment  waves  knocked  the  external  antennas  from 
their  plates  on  top  of  the  antenna  poles. 

3.  Pre-processing  Quality  Control 

To  mitigate  sources  of  error  as  best  as  possible  we  implemented  several  methods 
of  quality  control.  To  prevent  waves  from  knocking  antennas  from  drifters  after  the  first 
deployment,  the  external  antenna  was  glued  to  the  top  of  the  antenna  plate  using  epoxy. 

To  minimize  the  error  introduced  by  human  interaction  with  the  data,  efforts  were 
made  to  devise  a  way  to  turn  the  logging  devices  on  and  off  without  opening  the 
waterproof  boxes  while  the  drifters  were  carried  by  people.  These  efforts  were 
unsuccessful.  However,  it  was  recognized  that  isolating  the  external  antenna  from  the 
satellite  signal  would  prevent  the  loggers  from  recording  data  and  achieve  the  same 
result.  Caps  were  designed  to  cover  the  external  antenna  that  successfully  isolated  the 
signal.  However,  the  caps  proved  cumbersome  in  the  surf  zone,  and  were  rarefy  used. 
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For  each  deployment,  boundaries  were  established  outside  of  which  we  would 
retrieve  the  drifters  and  discount  data.  If  drifters  were  retrieved  inside  of  those 
boundaries  we  made  every  effort  to  move  quickly  up  the  beach  past  high  water  to  make  it 
obvious  that  motion  on  the  remainder  of  that  track  was  not  the  result  of  currents.  When 
drifters  exited  the  surf  zone  they  were  retrieved  by  boat.  Since  the  boat  would 
occasionally  move  at  slow  speeds  that  could  be  mistaken  for  currents,  drifter  tracks  were 
also  truncated  once  they  left  the  surf  zone. 

Ultimately,  a  note-taker  proved  to  key  for  quality  control.  The  note-taker  was 
responsible  for  recording:  time  each  logger  was  started  and  stopped;  time  each  drifter 
entered  the  water,  was  removed  from  the  water,  or  exceeded  observation  boundaries;  time 
each  drifter  became  tangled  or  untangled  from  kelp;  and  (during  the  first  deployment) 
time  each  drifter  antenna  was  knocked  off  by  a  wave  and  replaced  by  swimmers.  During 
deployments  the  note-taker  had  no  other  responsibilities  due  to  the  importance  of  the  job. 
Since  it  would  be  difficult  for  one  person  to  keep  track  of  so  many  drifters,  those 
deploying  drifters  helped  the  note-taker  by  reporting  observations  of  any  of  the  listed 
occurrences  to  the  note-taker. 

B.  DATA  PROCESSING 

1.  Post-processing 

The  DGPS  data  collected  by  the  drifters  were  uploaded  from  drifter  and  base 
station  logging  devices  and  post-processed  using  commercial  software.  Since  the  base 
station  was  stationary,  an  estimate  of  the  pseudo-range  bias  (difference  between  actual 
and  measured  distance  between  GPS  satellite  and  receiver,  based  on  upper-atmosphere 
effects)  was  calculated  for  each  data  point  from  the  base  station  (MacMahan  and  Brown, 
2007).  Since  the  base  station  and  drifter  antennas  were  in  close  proximity,  the  same 
upper-atmospheric  conditions  affect  them,  and  the  base  station  bias  for  each  data  point 
was  removed  from  the  drifter-logged  for  each  corresponding  time.  Position  errors  using 
non-differential  GPS  are  on  the  order  of  10m.  Use  of  DGPS  reduces  that  error  to  the 
order  of  centimeters  (MacMahan  and  Brown,  2007).  Processed  logger  files  were 
converted  to  text  format. 
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Drifter  positions  were  recorded  every  two  seconds.  In  order  to  make  the  positions 
more  meaningful,  the  global  Universal  Transverse  Mercator  (UTM)  coordinates  for  each 
observation  were  first  transformed  to  local  RCEX  coordinates.  Corresponding  u-  and  v- 
components  of  velocity  are  calculated  by  a  forward  differencing  scheme. 

2.  Post-processing  Quality  Control 
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Figure  13.  Unfiltered  time-series  data  for  a  single  drifter  on  5  May  07:  (a)  unfiltered 
cross-shore  position;  (b)  unfiltered  cross-shore  velocity.  Time  series  are  mostly 
smooth,  but  noise  in  the  positional  time-series  results  in  a  noisy  velocity  time- 
series. 


Once  the  data  were  in  local  coordinates  and  time,  deployment  notes  were  used  to 
remove  human-  and  kelp-influenced  data.  The  resulting  time-series  (Figure  13)  were 
overall  smooth,  but  still  had  obviously  erroneous  spikes  in  position  and  hence  velocity. 
To  remove  the  spikes,  the  data  were  filtered  by  removing  outliers  that  were  greater  than 
three  standard  deviations  from  the  mean  velocity  for  that  drifter  deployment.  After 
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filtering  twice  data  were  visually  checked  and  points  were  individually  removed  if  they 
met  velocity  criteria  but  still  exceeded  reasonable  spatial  separation. 


The  resulting  gaps  in  position  data  are  filled  using  a  spline  interpolation  for  gaps 
of  10s  or  less,  and  linear  interpolation  for  larger  gaps.  The  resulting  data  are  much 
smoother  than  before  the  quality  control  process.  (Figure  14) 
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Figure  14.  Time-series  data  after  quality  control,  for  the  same  drifter  as  Figure  13:  (a) 
cross-shore  position;  (b)  cross-shore  velocity.  Smoothness  is  significantly 
improved  after  filtering. 


3.  Velocity  Binning 

The  u-  and  v-components  of  velocity  are  grouped  by  dividing  the  field  site  into  a 
spatial  grids  on  5mx5m  and  10/77x10/77  bins  for  comparison.  The  means  and  standard 
deviations  of  the  velocity  components  are  calculated  for  each  bin.  Velocity  averages  are 
then  plotted  as  vectors  against  a  background  of  bathymetric  contours  to  show  a 
qualitative  relationship  between  the  flow  pattern  and  bathymetry  (Figure  1 5). 
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Figure  15.  Bin-averaged  velocity  vectors  plotted  against  bathymetry  using:  (a) 

5m  x  5m  grid  spacing;  and  (b)  10 m  x  10 m  grid  spacing.  The  average  shoreline  for 
the  deployment  period  is  represented  by  the  border  between  yellow  contours 
(sand)  and  blue  contours  (water).  Arrows  are  vectors  of  averaged  velocities 
within  each  bin.  Vectors  shoreward  of  the  shoreline  are  the  result  of  tidal  changes 
in  the  mean-water  level  during  the  deployment  and  possibly  from  imperfect 
removal  of  points  during  quality  control.  The  smaller  bin  size  provides  finer 
resolution,  but  introduces  more  noise,  while  the  larger  bins  provide  smoother 
profiles  at  the  expense  of  slightly  lower  resolution. 
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4.  Bathymetric  Adjustment 


While  underlying  beach  morphology  at  the  field  site  changes  little  over  time,  tides 
still  change  the  bathymetric  mean  sea  level  (0 m  contour)  with  time.  To  obtain 
bathymetry  relevant  to  a  particular  deployment  period,  an  adjustment  is  made  to 
bathymetric  surveys.  Hourly  mean  sea  levels  for  each  PUV  are  computed  relative  to  the 
same  0 m  contour  used  in  the  bathymetric  survey.  These  hourly  averages  are  further 
averaged  over  a  deployment  period  of  interest.  Bathymetry  for  the  most  relevant  survey 
is  adjusted  by  this  second  average  to  obtain  bathymetry  for  a  given  deployment. 
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IV.  RESULTS 


A.  RIP  CHARACTERISTICS 

To  determine  rip  current  and  rip  channel  characteristics,  cross-shore  velocity  and 
bathymetric  profiles  for  rip  channels  are  averaged  daily  and  for  the  entire  experiment. 
Alongshore  velocity  and  bathymetric  profiles  are  also  averaged  for  the  entire  experiment. 

1.  Rip  Currents 

Cross-shore  velocities  down  the  center  of  a  rip  channel  (Figure  16)  increase  to  a 
maximum  of  50 cm/s  at  about  50 m  from  the  shore.  Beyond  50m,  the  velocities  decay  to 
about  20 cm/s,  with  another  small  peak  at  about  100m,  after  which  they  decay  again.  The 
peak  in  velocities  at  100m  is  likely  a  result  of  fewer  observations  and  a  bias  towards 
higher  velocities  as  drifters  get  caught  in  rip  pulsations  that  eject  them  from  the  surf  zone. 
Average  rip  current  velocity  is  3 0cm A. 


Average  cross-shore  profiles  of  U  in  rip  currents.  Grid:  10m 


Figure  16.  Cross-shore  velocities  in  a  rip  channel.  Cross-shore  axis  is  normalized  for 
the  experiment  to  represent  meters  from  the  shoreline  rather  than  meters  from  the 
local  coordinate  system  alongshore  axis. 
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2.  Rip  Channels 


Cross-shore  bathymetric  profiles  in  the  center  of  rip  channels  show  increase  in 
depth  to  about  1.5 m  at  50 m  from  the  shoreline.  The  slope  for  the  first  50 m  is  about  1:33. 
The  profile  from  50 m  to  100/77  is  roughly  flat,  before  becoming  steep-sloped  again  after 
100/77  (Figure  17). 


Cross-shore  profile  of  rip  channel 


Figure  17.  Cross-shore  bathymetric  profile  in  a  rip  channel.  Cross-shore  axis  is 

normalized  for  the  experiment  to  represent  meters  from  the  shoreline  rather  than 
meters  from  the  local  coordinate  system  alongshore  axis. 


Comparison  of  rip  channel  cross-shore  velocity  profiles  to  rip  channel 
bathymetric  profiles  reveals  a  correlation  between  the  maximum  velocity  and  the  change 
in  slope  at  50 m.  The  increase  of  slope  at  100/77  also  corresponds  well  to  the  edge  of  the 
surf-zone  (e.g.,  breaker  location)  as  found  in  time-lapsed  imagery.  With  the  apparent 
relationship  between  surf  zone  width  and  increasing  slope  leads,  the  end  of  the  rip 
channel  is  defined  as  the  end  of  the  surf-zone.  From  this  definition,  the  average  channel 
length  is  100m,  and  the  average  maximum  channel  depth  is  1.7m. 
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Figure  18.  Alongshore  bathymetric  profile  ~5 Om  from  the  shoreline,  showing  how 
rip-channel  spacing  and  width  were  determined. 

Alongshore  bathymetric  profiles  are  used  to  determine  average  width  and  spacing 
of  rip  channels  (Figure  18).  The  profiles  used  are  about  50/77  from  the  shoreline,  where 
the  maximum  rip  current  velocities  are  found.  Rip  channel  width  is  measured  as  the 
distance  between  greatest  change  in  slope  between  the  channel  itself  and  the  shoal  on 
either  side  of  the  channel  -  the  average  width  is  55 m.  The  channel  spacing  is  the  distance 
between  depth  maxima  of  adjacent  rip  channels  -  average  channel  spacing  is  about  90 ///. 

B.  VOLUME  BALANCE 

Drifter  data  for  4  May  shows  gyres  in  surf-zone  circulation  centered  at  95 m  and 
105/77  offshore  in  local  coordinates.  Since  10m  bins  provide  a  smoother  velocity  profile 
than  5/7/  bins,  the  larger  bins  were  used  for  calculations  of  volumetric  flow  rate  (Figure 
19).  The  bin  velocities  were  further  smoothed  by  averaging  over  four  bins  (from  80/h  to 
1 10/// )  in  the  cross-shore  direction. 
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Figure  19.  Alongshore  profile  of  cross-shore  velocities  -50 m  from  the  shoreline. 

Positive  velocities  are  off-shore,  while  negative  velocities  are  onshore. 

To  verify  continuity,  a  bathymetric  profile  is  taken  along  a  transect  -50 m  from  the 
shoreline.  Discrete  volumetric  flow  (V)  is  then  calculated  along  the  transect  by: 

V  =  uh ■  Ay 

where  u  is  the  bin  average  cross-shore  component  of  velocity  for  the  given  point,  h  is 
depth  of  the  water  column  at  that  point,  and  Ay  is  the  bin  width.  This  assumes  no 
vertical  shear  in  the  water  column.  Summing  the  discrete  volumetric  flow  rate  over  two 
gyres  (from  30m  to  220m  alongshore)  gives  an  offshore  discharge  of  -58.62 m/s  and 
onshore  discharge  of  -58.56 rn/s.  Net  discharge  is  ~0.0646m5A,  which  is  an  error  of  less 
than  0.1%  considering  the  gross  discharge  is  ~117m  /s  (Figure  20).  Performing  similar 
balances  for  profiles  on  5  May,  7  May,  and  10  May  yields  an  overall  average  net- 
discharge  imbalance  of  about  6.67%  with  a  standard  deviation  of  4.79%. 


28 


Figure  20.  Volumetric  flow  rate  profiles  along  transect  located  ~50 m  from  the 
shoreline  -  positive  is  offshore  flow. 

The  peaks  in  Figure  20  roughly  correspond  to  rip-channel  centers  while  valleys 
correspond  to  shoals.  As  expected  from  continuity,  the  flow  rate  over  the  shoals  is 
characterized  by  a  relatively  flat  valley  along  the  width  of  the  shoal,  while  flow  in  the  rip 
channel  is  characterized  by  a  narrower  peak  of  greater  magnitude. 

C.  EDDY  CHARACTERISTICS 

Drifter  data  for  each  day  is  sectioned  into  individual  deployments.  Each  of  those 
individual  deployments  are  plotted  and  examined  to  find  those  with  enough  oscillatory 
motion  to  facilitate  a  general  description  of  circulation  cell  characteristics  for  each  day 
and  for  the  site  overall.  The  time-series  of  cross-shore  drifter  positions  are  demeaned, 
detrended,  and  auto-correlated  for  each  deployment.  At  a  lag  of  r  =  Os  the  auto¬ 
correlation  function  of  the  demeaned,  detrended  data  is  equal  to  the  variance 
( R  |  =  oy ),  which  permits  calculation  of  the  standard  deviation.  Since  95%  of  a 

dataset  is  contained  within  two  standard  deviations  of  its  mean,  the  eddy  radius  is  defined 
as  2<tv  .  Since  diameter  is  twice  the  radius,  eddy  diameter  is  4 crx .  Period  of  rotation 

corresponds  to  the  first  lag  ( ^  Os  )  where  the  auto-correlation  function  is  a  maximum. 
(Figure  21;  Table  1) 

Based  on  analysis  of  drifter  data  from  days  on  which  eddies  were  observed,  the 
mean  eddy  diameter  for  the  experiment  is  88/u  and  the  mean  period  is  4.7 min.  On 

average,  each  drifter  experienced  at  least  eight  revolutions  in  the  surf  zone  before  landing 
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on  shore  or  heading  outside  the  surf  zone.  The  consistency  of  these  data  and  their 
correspondence  to  surf  zone  characteristics  (e.g.,  surf  zone  width)  indicate  that  these 
eddies  are  real  and  persistent. 


Figure  2 1 .  Sample  autocorrelation  plot  of  cross-shore  positions.  Lag  (seconds)  is  the 
x-axis  and  auto-correlation  value  (nr)  is  the  y-axis.  Since  the  auto-correlation  is 
based  on  demeaned  and  detrended  data,  the  peak  at  zero  lag  represents  the  square 
of  the  data  variance.  The  value  of  lag  at  the  next  peak  in  auto-correlation 
corresponds  to  the  period  of  oscillation  for  the  data  set. 
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Average 

Diameter 

Period 

27  Apr  07 

92.0  m 

267s  (4.5 min) 

4  May  07 

86.0m 

255 s  (4.3 min) 

5  May  07 

91.6  m 

285 s  (4.8 min) 

7  May  07 

82.8  m 

319s  (5.3 min) 

10  May  07 

16.9  m 

240 s  (4.0mm) 

15  May  07 

12.9  m 

307s  (5.1mm) 

Overall 

Average 

8  8 . 0m 

282s  (4.7mm) 

Table  1.  Eddy  characteristics  for  the  six  deployment  days  with  any  observed  eddies.  The 
“overall”  average  excludes  measurements  from  10  May  and  15  May  -  even 
though  eddies  were  observed  on  those  days,  flow  was  primarily  alongshore. 

D.  OBSERVATION  DENSITY 

An  interesting  feature  of  the  eddies  is  their  fairly  persistent  center  positions.  The 
highest  number  of  observations  occurs  near  the  center  of  each  circulation  pattern  on  an 
observation  density  plot  (Figure  22).  This  indicates  that,  during  the  deployment,  drifters 
spend  a  significant  amount  of  time  near  the  eddy  centers  -  observations  are  2s  apart  and 
average  velocities  are  lower  at  the  eddy  centers  than  the  edges.  This  is  an  interesting 
phenomenon  that  could  have  an  impact  on  pollution  or  biota  transport  (Talbot  and  Bate, 
1987). 
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Drifter  Survey  050507,  Grid:  10m 
Observation  Density 


alongshore,  m 

Figure  22.  Observation  density  plot,  showing  higher  concentration  of  drifter 
observations  at  the  centers  of  surf  zone  eddies. 

E.  EULERIAN  COMPARISON 

To  validate  drifter  observations  they  are  compared  with  in  situ  instruments  that 
are  already  widely  accepted  as  correct.  Hourly  means  of  the  PUY  and  ADCP  recorded 
speeds  are  compared  with  corresponding  hourly  means  of  drifter  speeds  within  10m  bins 
centered  on  each  instrument  in  the  array.  Bins  with  fewer  than  30  observations  are 
filtered  out  since  they  are  statistically  insignificant.  A  scatter  plot  of  the  remaining 
speeds  reveals  a  “shot-gun”  pattern;  however  the  linear  regression  is  close  to  a  one-to-one 
fit  with  an  r2  of  0.81  (Figure  23). 
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Hourly-avg  speed  of  drifters  Win  5m  of  PUVs 
vs  Hourly-avg  PUV  speed  For  all  days. 


Figure  23.  Results  of  a  linear  regression  comparing  average  in  situ  instrument 

velocities  with  average  drifter  velocities  near  them.  While  there  is  scatter  in  the 
data  correlation  is  relatively  high  (r-squared  of  0.81,)  with  a  best  fit  slope  just  less 
than  one,  and  y-intercept  only  slightly  larger  than  zero. 
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V.  CONCLUSIONS 


A.  SUMMARY 

1.  General  Characteristics 

Rip  current  circulation  cells  are  found  to  be  contained  within  the  surf  zone.  The 
surf  zone  is  the  energetic  region  that  extends  from  the  shoreline  to  the  offshore  location 
where  waves  begin  to  break.  The  beach  profile  transitions  from  a  low-tide  terrace  (slope 
-1:100)  to  an  offshore  slope  (about  1:20)  approximately  100m  from  the  shoreline,  which 
is  the  surf  zone  width.  Rip  channel  length  scales  are  characterized  by  an  alongshore 
spacing  of  O(90m),  cross-shore  channel  length  of  0(1 00m),  and  channel  width  of 
0(55m).  The  morphodynamic  relationship  for  this  site  is  governed  by  wave  breaking  and 
bathymetry,  generating  patterns  of  onshore  flow  over  shoals  and  offshore  flow  in  narrow 
channels,  rip  currents. 

Cross-shore  velocities  in  the  rip  current  increase  from  zero  at  the  shoreline,  to  a 
maximum  of  O(60cm/s  )  midway  through  the  surf  zone,  and  decay  to  the  edge  of  the  surf 
zone.  An  additional  smaller  peak  in  velocities  was  observed  to  occur  at  the  edge  of  the 
surf  zone,  and  is  believed  to  be  related  to  rip  current  pulsations.  Average  rip  current 
velocities  are  30  ±  15cm  /  s  for  seven  deployments  at  various  wave  and  tidal  conditions. 

Drifters  deployed  during  RCEX  demonstrate  the  presence  of  persistent  surf  zone 
eddies  characterized  by  semi-closed  loops  related  to  the  rip  current  circulation.  These 
eddies  complete  a  rotation  in  about  4.7mm,  and  are  about  88m  in  diameter.  Continuity  of 
shoreward  and  seaward  discharge  calculations  about  50m  from  shore  (near  the  eddy 
centers)  balance  to  less  than  10%  of  the  gross  discharge. 

2.  Drifter  Performance 

With  ground-breaking  success,  a  low-cost  drifter  has  demonstrated  its  ability  to 
significantly  improve  our  picture  and  understanding  of  the  horizontal  flow  structure  near 
shore  and  fill  in  gaps  left  by  using  only  in  situ  instruments.  Comparison  of  average 
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drifter  speeds  near  in  situ  instruments  shows  relatively  high  correlation  between  the  two 
measurements.  The  r-squared  value  of  -0.81  for  a  linear  regression  fit  between  the  two 
sets  of  data  demonstrates  that  drifter  speeds  compare  well  with  speed  obtained  by  in  situ 
instruments  that  are  already  accepted  as  accurate. 

3.  Setting  the  Record  Straight 

The  prevailing  opinion  is  that  the  best  way  out  of  a  rip  current  is  to  swim  parallel 
to  shore  -  this  is  not  necessarily  good  advice.  The  presence  of  eddies  on  either  side  of  a 
rip  current  suggest  that  a  swimmer  would  at  best  have  a  50%  chance  of  swimming  with 
the  current  by  swimming  parallel  to  the  shore.  A  swimmer  has  a  better  chance  of 
conserving  energy  and  returning  to  land  safely  by  floating  with  the  current  and  letting  it 
carry  him  back  to  shore  with  the  onshore  flow  of  the  eddy.  In  doing  so,  the  swimmer  will 
still  need  to  deal  with  oncoming  waves,  especially  as  he  nears  the  edge  of  the  surf  zone, 
but  he  will  expend  less  energy  in  the  process. 

B.  OTHER  WORK 

1.  Current  Research 

Since  RCEX  only  took  place  six  months  ago,  this  thesis  contains  the  preliminary 
results  of  a  much  larger  picture  of  rip  current  circulation  dynamics.  Other  research  in 
progress  based  on  the  experiment  includes:  measurement  of  dispersion  and  diffusivity 
based  on  multi-particle  statistics  through  the  drifters;  investigation  of  the  effect  of 
waves/tides  on  circulation  cells  and  possible  rip  current  pulsations;  and  examination  of 
the  poor  correlation  between  individual  velocity  components  of  the  Lagrangian  and 
Eulerian  measurements. 

2.  Future  Research  and  Improvement 

Other  interesting  research  that  could  stem  from  this  are:  exploration  of  the  role 
these  circulation  patterns  play  in  sediment  transport;  or  examination  of  the  effect  of  these 
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nearshore  circulation  patterns  on  humans  by  making  similar  measurements  using  human 
“drifters.”  It  would  also  be  useful  to  extend  these  results  to  modeling  in  order  to  better 
predict  flow  behavior  near  shore. 

Future  experiments  similar  to  RCEX  could  be  improved  by  finding  some  way  to 
turn  the  logging  devices  on  and  off  without  having  to  open  the  water-proof  boxes  during 
deployment,  thus  reducing  error  introduced  by  humans.  It  might  also  prove  useful  to 
have  a  compass  as  an  integral  part  of  the  in  situ  instruments  to  eliminate  human  error 
associated  with  determining  their  orientation. 

3.  Application 

The  results  presented  here  have  a  direct  application  to  the  Navy.  SEAL  and  mine¬ 
hunting  operations  both  deal  directly  with  the  surf-zone.  An  understanding  of  the  rip 
current  circulation  patterns  in  the  surf  zone  will  improve  predictive  modeling  and  aid  in 
planning  mine-hunting  and  SEAL  delivery. 
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